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A new flare star member candidate in the Pleiades cluster 
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We present a new flare star, which was discovered during our survey on a selected field at the edge of the Pleiades cluster. 
The field was observed in the period 2007 - 2010 with three different CCD-cameras at the University Observatory Jena 
with telescopes from 25 to 90 cm. The flare duration is almost one hour with an amplitude in the R-band of about 1.08 
mag. The location of the flare star in a color-magnitude diagram is consistent with age and distance of the Pleiades. In the 
optical PSF of the flare star there are two 2MASS objects (unresolved in most images in the optical Jena PSF), so it is not 
yet known which one of them is responsible for this flare. The BVRIJHK colors yield spectral types of Ml and M2 with 
extinction being Av = 0.23 1 ± 0.024 mag and Ay = 0.266 ± 0.020 for those two stars, consistent with the Pleiades cluster. 
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1 Introduction 

The first flare star was discovered in 1924 by Hertzsprung 
(Hertzsprung, 1924). Later a new type of flares was dis- 
covered in the solar neighborhood and carried the name 
UV Cet-stars, which are known to be low-mass stars with 
very low luminosities. The investigation and study of the 
flares in nearby open clusters such as the Pleiades showed 
their importance for the evolution of young low-mass stars 
and protostars (Mirzoyan et al. 1995). The first flare star in 
the Pleiades was discovered by Johnson & Mitchel (1958). 
By studying the statistical characteristics of flare stars in 
the Pleiades by Ambartsumyan et al. (1970, 1971), it was 
found that almost all stars fainter than V = 13.3 mag can be 
flare stars. Another advantage of flare star investigations in 
nearby clusters is to identify additional cluster members at 
very low masses. Most of the Pleiades faint members were 
discovered in flare star surveys (Jones 1981; Haro et al. 
1982; Stauffer et al. 1991). 

A field at the edge of the Pleiades cluster was chosen to 
be observed in our survey to find new low-mass stars or 
brown dwarfs and variable stars, because in this field there 
was no survey for brown dwarfs done before (see Table 1 
in Schwarz & Becklin 2005). One new flare star was found, 
which lies on the Pleiades zero-age main-sequence, and is 
presented in this paper. In Sect. 2 we describe our observa- 



* Corresponding author: e-mail: mmoualla@astro.uni-jena.de 



tions and data reduction. The photometric analysis of our 
data is shown in Sec. 3, results are presented and discussed 
in Sect. 4. 



2 Observations and data reduction 

Our field was chosen to be at the edge of the Pleiades cluster, 
in order to look for new brown dwarfs and low-mass stars 
by deep imaging and follow-up spectroscopy, since no deep 
CCD surveys were done there before; at the same time, we 
wanted to study the variability of all stars in the field on 
timescales of minutes to years. Concerning brown dwarfs 
several candidates were found in our field (Eisenbeiss et al. 
2009; Seeliger et al. 201 1, in preparation). 
We observed two overlapping sub-fields with field centers 
at ajsooo = 3 h 42 m 20.6 s , £j 20 oo = +25° 36' 54" and also 
at a J200 o = 3 h 40 m 54 s , <5 J200 o = +25° 40' 32". 
All observations were carried out at the University Observa- 
tory Jena, which is located close to the village GroBschwab- 
hausen (GSH), about 10 km to the west of Jena. In March 
2007 we started our survey using the 0.25 m Cassegrain- 
Teleskop-Kamera (CTK), which is installed at the 0.25 m 
auxiliary Cassegrain telescope, mounted at the tube of the 
0.9 m telescope (see Mugrauer 2009). It has a field of view 
(FoV) of 37.7' x 37.7'. In August 2010 the CTK was re- 
placed with a new CCD-camera called CTK-II, which has 
a FoV of 21.0' x 20.4' (Mugrauer 2011, in preparation). 
Since beginning of February 2009 the new Schmidt-Teleskop-| 
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Kamera (STK) is in operation in the Schmidt-focus of the 
0.9 m telescope (0.6 m in Schmidt mode) at the University 
Observatory Jena. With its 2048x2048 pixels CCD detec- 
tor and a pixel-scale of about 1.5 arcsec/pixel it has a FoV 
of 52.8' x 52.8' (see Mugrauer & Berthold 2010). Since 
the STK camera was used mainly in our survey due to its 
larger FoV in comparison to the other cameras we list in 
Table |2]only the observations with the STK camera. Obser- 
vations were done mostly in R-band. At the beginning of 
our survey with the CTK we used an individual exposure 
time of 60s per image. The same exposure time was used 
later with the STK. We found many new variable stars in 
the field (Moualla 201 1). One of them was very faint, so that 
we increased our exposure time to 90 s. Basic data reduc- 
tion was done for all images, by means of dark subtraction 
and flat fielding using the ESO software MIDAS (Munich 
Image Data Analysis System). 
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Fig. 2 A zoom ( 1 ' x 1.2') into the FoV of the STK camera 
shows that the new flare star, observed in R-band with an 
exposure time of 90s, is one of two stars (unresolved in most 
images, so that individual light curves cannot be obtained). 
North is up, East is to the left. 



3 Photometry 

After reducing the images, three different programs were 
used to obtain positions and magnitudes of all objects in the 
field. 

First the source extractor (SE) for the source detection with 
GAIA (Graphical Astronomy and Image Analysis Tool) was 
performed on each single frame in order to distinguish be- 
tween real and false objects by performing thresholding and 
deblending (see Bertin 1997) as well as to determine their 
center coordinates (Xc cn , Yc en )- 

SE uses automatic aperture photometry, originally designed 
for galaxies and extended objects; especially in crowded 
fields as well as undersampled point sources, this method 
fails to obtain accurate magnitudes; furthermore, their er- 
rors are underestimated; however, source detection based on 
thresholding, as SE employs, is a fast and convenient way 
to detect even faint objects with reasonable positional accu- 
racy. 

The instrumental magnitudes for all objects detected by SE 
were determined then by using aperture photometry with 
MIDAS. We created a program for the aperture photome- 
try with MIDAS which needs the center coordinates of all 
objects in the field. The flux of the star was measured with 
an internal aperture with a radius of 5 pixels and the back- 
ground was measured 20 pixels away around the internal 
one. 

The differential magnitudes of all objects in the field were 
determined by using an algorithm (Broeg et al. 2005) that 
uses all stars in the field as comparison stars in order to 
determine the best synthetic comparison star by calculat- 
ing their weighted average. Very small weights were given 
to unsuitable stars that are possibly variable or very faint. 
By excluding these unsuitable stars the best synthetic com- 
parison star is calculated and the differential magnitudes of 
all objects in the field are determined by comparing their 
instrumental magnitudes to the magnitude of the synthetic 
comparison star. 



Fig. 3 A zoom (33.80" x 22.62") into the Ks-band 
2MASS image of our flare star (Skrutskie et al. 2006). North 
is up, East is to the left. Stars 1 and 2 as in Table 5. 

To determine the absolute magnitude of the field objects we 
used those well-known stars in our field that are constant 
(from their light curves in all observed nights). Those stars 
are listed in Table [T] The apparent magnitude m can be cal- 
culated by the following equation 

m = C + m instr - k • z, (1) 

i.e. for the R-band we get 

R = C R + Rinstr - k R • z (2) 

where R is the magnitude to be measured and Rinstr refers 
to the measured instrumental magnitude with our camera, 
k and z are the extinction coefficient and airmass, respec- 
tively, whereas Cr denotes the zero point of the detector 
in the R-band. Long-term variations of color terms are not 
yet available for this instrument, but should also be negligi- 
ble for short-term events like flares. For the determination 
of the zero point we used the constant known bright stars 
(Table Q]) in our Pleiades field, so that all stars in the field 
have the same extinction coefficient (k) due to the same air- 
mass (z). To calibrate our images we added the zero point 
value (mean value in Table 1) to the measured differential 
magnitude (m - C + k • z) of every single star in the field. 
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Fig. 1 The left panel shows the Pleiades field with the flare star observed by the STK camera in R-band in the night 2010 
September 20 with an exposure time of 90s. The field center on the left is a,/ 2 ooo = 3 h 40 m 54 s , <5,/ 2 ooo = +25° 40' 32". 
The right panel zooms into the FoV of the STK, the circle marks the new flare star. The FoV of the left image is 52.8' x 
52.8' and the right field is 3.9' x 2.9'. North is up, East is to the left. 

Table 1 List of constant stars in the field used to determine the zero point of the detector in VRJ bands. In fall 2010 the 
V-filter showed impurity, so that the quality of the measurements in V is lower than in R and I. 



Star name 


RA(J2000.0) 


Dec(J2000.0) 


(Cv-k-z) 


(Ca-k-z) 


(Ci-k-z) 


Spty 




[h m s] 


[° ' "] 




[mag] 






Melotte 22 SK 785 


03 40 14.80 


+25 19 20.0 


13.862±0.174 


13.533 ±0.022 


13.024 ±0.017 


M3V 


Melotte 22 MT 22 


03 40 58.40 


+25 17 43.0 


13.357 ±0.268 


13.297 ±0.016 


12.402 ±0.013 


A7V 


Melotte 22 SK 792 


03 40 05.00 


+25 31 38.0 


13.960 ±0.216 


13.083 ±0.019 


12.671 ±0.009 


K4V 


NVSS J034035+253859 


03 40 35.57 


+25 39 00.6 


13.347 ±0.180 


13.608 ±0.015 


12.365 ±0.018 


A8V 


Melotte 22 DH 153 


03 40 49.89 


+26 03 33.5 


13.935 ±0.191 


13.387 ±0.014 


12.626 ±0.009 


M2V 


Melotte 22 SK 769 


03 40 24.60 


+26 03 26.0 


13.302 ±0.262 


13.121 ±0.018 


12.591 ±0.014 


K4V 


Mean value 






13.627 ±0.215 


13.338 ±0.0171 


12.613 ±0.013 





4 Results 

Our flare star is one of two unresolved stars (Fig. 1 and 
2, unresolved in most images) with similar brightnesses in 
JHK-bands. In the 2MASS image (Fig. [3}, the two stars are 
well resolved. All information about both components are 
gathered in Table [5] The coordinates and the magnitudes 
in infrared are taken from the 2MASS catalog (Cutri et al. 
2003). The B magnitude (B = 18.270mag) is from the NO- 
MAD catalog (Zacharias et al. 2005). Since we are unable 
to resolve the two components, the measured magnitudes 
in VRI-bands belong to both components together (Table 
|4j. In order to calculate the individual magnitudes for each 
component we first used the individual JHK magnitudes of 
every component to determine its spectral type taking into 
account the optical extinction (Ay) of the Pleiades of 0.2 
mag (van Leeuwen 2009), which was used to calculate the 
extinction in other bands using Table 3 in Rieke and Lebof- 
sky (1985). Then, by using the spectral type, the JHK mag- 
nitudes and Table A5 in Kenyon and Hartmann (1995) we 
determined the BVRI magnitudes of each component indi- 
vidually, which can be used to calculate the flux ratio of 



both components in BVRI-bands. Finally, we determine the 
individual magnitudes of each component using the flux ra- 
tio and the measured magnitudes of both components to- 
gether. The subtracted extinction is added again to all mag- 
nitudes, because we calculate the final spectral type and ex- 
tinction by using the BVRIJHK magnitudes and Table 3 in 
Rieke and Lebofsky (1985) and Table A5 in Kenyon and 
Hartmann (1995). The final spectral type of both stars (Ta- 
ble 5) is Ml-2 and the derived extinction is in agreement 
with that of the Pleiades. The effective temperature can be 
easily determined using Table A5 in Kenyon and Hartmann 
(1995) since the spectral type is known. Mass and radius are 
taken from Schmidt-Kaler (1982) assuming that both stars 
are main sequence stars. 

On the night 2010 September 22, a flare was seen in the 
combined light curve of these two stars, whereas they show 
no variability but only a constant signal in all other ob- 
served nights. The light curve is typical for flares: Constant 
light before and after the flare, a short decrease in brightness 
(negative flare) before the sharp increase, then an exponen- 
tial decay. 

In Fig. [5] we present the observed light curve of the flare, 
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Fig. 4 



Light curves of the constant stars in the field used to determine the zero point of the detector in R-band. They 
were observed in the night 2010 October 21 with the STK in R-band using an exposure time of 90s. Shown from top to 
bottom are Melotte 22 SK 785 with a = 0.0163 mag, Mellote 22 MT 22 with a = 0.0057 mag, Mellote 22 SK 792 with a 
= 0.0165 mag, NVSS J034035+253859 with a = 0.0073 mag, Mellote 22 DH 153 with a = 0.0065 mag and Melotte 22 SK 
769 with cr = 0.0064 mag (see TableQ}. 



and its parameters indicated in Fig.[5]are given as values in 
Table [3] We have first fitted relative fluxes (flare flux rela- 
tive to the comparison star) of the flare decay part using the 
formula 



Flux = Const + A ■ exp(B ■ t) 



(3) 



The flux, the constant, and the factor A are unitless, because 
we fit a relative flux, i.e. the flux of the flare star relative to 
the synthetic comparison star; time t is in days, constant B 
has then the unit 1/day. With this non-linear model we ob- 
tained a very good fit with parameters: 
A = 0.0637 ±0.0022, 
B = -189.64 ±10.86 1/day, 
Const = 0.04219 ±0.00031. 

In order to determine the flare end time, in addition to the 
above mentioned model, we considered also purely constant 
as an alternative one. Thus, we compared two models via F- 
statistics, i.e. the null hypothesis that the first model (expo- 
nential decay of the flare) does not provide a significantly 
better fit than second one (purely constant). We required 
the false rejection probability of the null hypothesis to lie 



below 0.05, the usual value. We applied this approach sub- 
sequently to all data subsets starting from the flare maxi- 
mum up to the end of the observations. The purely constant 
model can be accepted (rejecting the flare model) from the 
time JD-2455461 = 0.5913 to the end of observations, i.e. 
the flare decay time lasted 1 hour, which is typical for flares 
in the Pleiades cluster; from that time on, the probability for 
the emission to be constant again (after the flare) is always 
> 0.95. 

Within the total exposure time of 132 hours we observed 
one flare with an amplitude of about 1 .08 mag in the R-band. 
Haro et al. (1982) have observed in the Pleiades a total of 
1800 flares (in the blue wavelength range with an average 
amplitude of about 2.5 mag) in 519 stars during a total ex- 
posure time of 3500 hours, i.e. one flare on average every 
972 hours. 

We have observed in our field 11-12 Pleiades flare stars (11 
known and one probably new flare star of the Pleiades). 
Thus on average every (972 / 12 =) 81 hours a flare is ex- 
pected - for a total exposure time of 132 hours thus 1.6 
flares. This is not inconsistent with the actual observed flare. 
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Fig. 5 Light curve with the flare observed in the night 2010 September 22 with the STK camera in R-band with an 
exposure time of 90 s. We plot relative flux, i.e. flare flux relative to the synthetic comparison star, which has 13.338 mag. 
The flare began at 1:07(3) UT (JD 2455461.536), reached its maximum with an amplitude > 1.08 mag at 01:11(1.5) UT 
(JD 2455461.544), and ended at 02:03(4)UT (JD 2455461.581) with a total duration of 56 ±4min. (see Table [5j. The 
full line is the best fit to the flare decay, see text. The significance for the detection of a short small negative flare at JD 
2455461.545 is only 2.6 a. 



None of the two stars (one of which has flared) is mentioned 
in the General Catalog of Variable Stars (GCVS) (Samus et 
al. 2009) as flare stars. 

The proper motion of the second star 2MASS J03400685+ 
2514256 is \i a ■ cos S = 43.9 ± 3.8 mas • (yr) -1 and fig = 
- 44.3 ± 3.8 mas • (yr)" 1 (PPMXL catalog (Roeser et al. 
2010)), which is inconsistent with the mean of that of the 
Pleiades members n a ■ cosS = 20.10 ± 0.28 mas • (yr) _1 
and fig = -45.39 ± 0.27 mas • (yr) _1 (van Leeuwen 2009 ). 
Given the location and proper motion of this star (second 
star) with respect to the Pleiades, it also cannot be an ejected 
Pleiades member. 

For the first star, 2MASS J03400665+25 14276, there are no 
proper motion measurements available (NOMAD catalog, 
Zacharias et al. 2005). However, it might be that the two 
stars were also unresolved in the observations used to de- 
termine the proper motion of the second star, so that this 
proper motion may not be reliable. 



The positions of the possible flare stars in the color magni- 
tude diagram (Fig.|6]and|7]) are consistent with many known 
Pleiades members in our field as well as with low-mass 
Pleiades members, which are presented with three of the 
least massive stars (HHJ 2, 3, 6) in the Pleiades cluster (Ham-( 
bly et al. 1993). They are separated from Pleiades BDs like 
Calar 3 and Teide 1 (Rebolo et al. 1996) by PP1 15 (Stauf- 
fer et al. 1994) which defines with its mass M = 0.08 M 
(Basri et al. 1996) the borderline between low-mass stars 
and BDs in the Pleiades cluster. 

Both stars are consistent with the age, the extinction and the 
distance of the Pleiades members. 

We note that the Siess et al. (2000) isochrones in Fig. 6 (op- 
tical) seem to lie somewhat too low compared to our stars, 
very low-mass stars and brown dwarfs plotted (with litera- 
ture data), while the Siess et al. (2000) isochrones in Fig. 7 
(infrared) seem to lie too high, this is probably due to the 
models for very low masses. 
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Fig. 6 The color-magnitude diagram (corrected for extinction) for The two data points well below the main sequence are 
probably due to variability and non-simultaneous observations in R and I. Pleiades members in our field (diamonds) and 
the two 2MASS stars (asterisks) in our field, one of which has shown the flare. Five theoretical isochrones represent five 
ages (10, 25, 50, 80 and 150 Myr) and are taken from Siess et al. (2000) after adjusting their magnitudes to the distance 
of the Pleiades stars. Two data points well below the main sequence are probably due to variability and non-simultaneous 
observations in R and I. The point in the lower left shows the mean error bar for our data for the Pleiades members. 



In order to confirm or reject the first 2MASS star, whose 
proper motion is not yet known, as true Pleiades member, 
we plan to measure its proper motion with a new epoch IR 
image with high spatial resolution and to take a spectrum 
to confirm its spectral type. Since it is an early M-type flare 
star, it is probably a Pleiades member. 
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Table 2 

STK. 



Observations log of the Pleiades field with the 



Date 




Interval 
UT 




Filter 
Bessel 


Expo, 
[s] 


N 

images 


2009 Sep 18 


01: 


:50- 


-03: 


:43 


R 


60 


96 


2009 Sep 19 


00: 


16- 


-03 


:45 


R 


60 


175 


2009 Sep 21 


22: 


:41- 


-01 


:51 


R 


60 


161 


2009 Sep 24 


01: 


:01- 


-03 


:51 


R 


60 


142 


2009 Sep 26 


02: 


:56- 


-03 


:59 


R 


60 


53 


2009 Okt 19 


20: 


18- 


-22 


:30 


R 


60 


111 


2009 Okt 20 


20: 


:05- 


-23 


:30 


R 


60 


178 


2009 Okt 30 


23: 


09- 


-04 


:50 


R 


60 


243 


2009 Nov 06 


00: 


:45- 


-02 


:55 


R 


60 


95 


2009 Nov 14 


01: 


:02- 


-04 


:33 


R 


60 


91 


2010 Feb 03 


20: 


:42- 


-22 


:12 


R 


60 


74 


2010 Feb 16 


18: 


:41- 


-22 


:53 


R 


60 


206 


2010 Feb 20 


20: 


:05- 


-22 


:30 


R 


60 


120 


2010 Feb 24 


19: 


:27- 


-22 


:13 


R 


60 


140 


2010 Mar 02 


22: 


:21- 


-23 


:09 


R 


60 


41 


2010 Mar 03 


20: 


:37- 


-22 


:00 


R 


60 


60 


2010 Mar 06 


20: 


:40- 


-20 


:52 


R 


60 


11 


2010 Mar 07 


18: 


13- 


-19 


:34 


R 


60 


58 


2010 Mar 23 


18: 


:44- 


-20: 


:49 


R 


60 


74 


2010 Aug 21 


01: 


:20- 


-03: 


:10 


R 


90 


72 


2010 Sep 11 


00: 


:54- 


-02 


:32 


V/I 


180/50 


10/10 


2010 Sep 13 


23: 


:21- 


-02 


:29 


R 


90 


111 


2010 Sep 15 


23: 


00- 


-02 


:31 


R 


90 


75 


2010 Sep 17 


22: 


:21- 


-03 


:55 


R 


90 


198 


2010 Sep 18 


00: 


:32- 


-03 


:50 


R 


90 


75 


2010 Sep 20 


21: 


:59- 


-04 


:11 


R 


90 


221 


2010 Sep 21 


21: 


:58- 


-04 


:11 


R 


90 


219 


2010 Sep 22 


21: 


:59- 


-03 


:53 


R 


90 


210 


2010 Sep 23 


01: 


06- 


-03 


:47 


R 


90 


96 


2010 Okt 05 


22: 


:09- 


-23 


:39 


R 


90 


54 


2010 Okt 08 


21: 


12- 


-01 


:00 


R 


90 


134 


2010 Okt 09 


01: 


:55- 


-03 


:21 


V/I 


180/50 


20/20 


2010 Okt 10 


20: 


:05- 


-22 


:40 


R 


90 


88 


2010 Okt 11 


21: 


06- 


-21 


:52 


R 


90 


28 


2010 Okt 13 


20: 


06- 


-22 


:41 


R 


90 


92 


2010 Okt 17 


20: 


:31- 


-03 


:14 


R 


90 


239 


2010 Okt 20 


20: 


:39- 


-01 


:12 


R 


90 


114 


2010 Okt 21 


20: 


:31- 


-04 


:24 


R 


90 


280 


2010 Okt 22 


20: 


:29- 


-02 


:57 


R 


90 


230 


2010 Okt 28 


19: 


:21- 


-03 


:02 


R 


90 


222 


2010 Nov 02 


20: 


00- 


-22 


:34 


R 


90 


79 


2010 Nov 09 


19: 


:05- 


-21: 


:37 


R 


90 


88 



N is the number of images. 

Table 3 Parameters of the flare derived from its light 
curve in the night 2010 September 22 (JD 2455461). 



Parameter 


Value 


JD-2455461 


negative flare [h:m] 


01:05(2) UT 


0.545 


beginning of flare [h:m] 


01:08(3) UT 


0.5475 


duration negative flare [min] 


8 ±2 




maximum of flare [h:m] 


01:12(1.5) UT 


0.550 


limit of increase [min] 


< 4 




duration of decrease [h:m] 


00:59(2.5) 




end of flare [h:m] 


02:11(2) UT 


0.5913 


total time [min] 


63 ±4 




amplitude [mag] 


>1.08 





Table 4 Measured magnitudes of the flare star by using 
our STK camera (both components together unresolved), 
2010 October 21. 



V 


R 


I 


[mag] 


[mag] 


[mag] 


17.302 ±0.290 


16.968 ±0.084 


14.844 ±0.041 



Table 5 Characteristics of our flare star: VRI from our GSH data (split 
correctly among the two stars, see text), the B magnitude is taken from 
the NOMAD catalog (Zacharias et al. 2005), position, name, J, H, and 
Ks from 2MASS, spectral types and extinction from BVRIJHK (see text), 
temperature, mass, and radius for main sequence stars of those spectral 
types. 



Parameter 


star 1 


star 2 


a [h m s] 


03 40 06.657 


03 40 06.857 


8 [° ' "] 


+25 14 27.65 


+25 14 25.63 


2MASS-name 


03400665+2514276 


03400685+2514256 


B [mag] 


18.766 ±0.032 


19.358 ±0.063 


V [mag] 


17.817 ±0.291 


18.359 ±0.296 


R [mag] 


17.533 ±0.089 


17.946 ±0.105 


I [mag] 


15.523 ±0.052 


15.674 ±0.075 


J [mag] 


14.368 ±0.072 


14.437 ±0.032 


H [mag] 


13.764 ±0.088 


13.833 ±0.050 


K [mag] 


13.536 ±0.074 


13.563 ±0.038 


Sptype 


M1V (M0...M2) 


M2V (M1...M4) 


A v [mag] 


0.231 ±0.024 


0.266 ± 0.020 


T off [K] 


3720 (3850...3580) 


3580 (3720...3370) 


M[M ] 


0.45 (0.51. ..0.40) 


0.40 (0.51...0.27) 


R [R©] 


0.55 (0.60...0.50) 


0.50 (0.55...0.34) 
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